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ABSTRACT 

'in 

^ ! We investigate heating of the cool core of a galaxy cluster through the dissipa- 

tion of sound waves excited by the activities of the central active galactic nucleus 
(AGN). Using a weak shock theory, we show that this heating mechanism alone 
^ cannot reproduce observed temperature and density profiles of a cluster, because 

the dissipation length of the waves is much smaller than the size of the core and 
thus the wave energy is not distributed to the whole core. However, we find that 
if it is combined with thermal conduction from the hot outer layer of the cluster, 
the wave heating can reproduce the observational results. 
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1. Introduction 

o 

H 

The radiative cooling time of gas in the central regions of galaxy clusters (cool cores) 
is generally much smaller than the Hubble time. In the absence of any heating sources, 
this means that the intracluster medium (ICM) flows subsonically toward the cluster center 
with a mass deposition rate of M ~ 100-1000 M yr _1 (Fabian 1994). This flow was called 
a "cooling flow". However, recent X-ray observations have shown that the cooling rate of 
the ICM is much smaller. The Japanese ASCA team indicated that metal emission lines 
from the low temperature cooling gas were much weaker than that had been predicted by 
the classical cooling flow model (Ikebe et al. 1997; Makishima et al. 2001). This has been 
confirmed by XMM-Newton (Peterson et al. 2001; Kaastra et al. 2001; Tamura et al. 2001); 
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the actual mass deposition rates are about 1/10 of those predicted by the classical cooling 
flow model. 

The lack of metal emission lines indicates that the gas is prevented from cooling by some 
heating sources. At present, the most popular candidate for the heating source is the active 
galactic nucleus (AGN) at the cluster center (Tucker & Rosner 1983; Rephaeli 1987; Binney 
& Tabor 1995; Ciotti & Ostriker 2001; Bohringer et al. 2002; Churazov et al. 2002; Reynolds, 
Heinz, & Begelman 2002). However, it is not understood how the energy ejected by the AGN 
is transfered into the surrounding ICM. One idea is that bubbles inflated by AGN jets move 
outward in a cluster by buoyancy and mix the surrounding ICM (Churazov et al. 2001; 
Quilis, Bower, & Balogh 2001; Saxton, Sutherland, & Bicknell 2001). As a result of the 
mixing, hot ICM in the outer region of the cluster is brought into, and subsequently heats, 
the cluster center. The other idea is that the dissipation of sound waves created through 
the AGN activities. In fact, sound waves or weak shocks that may have evolved from sound 
waves are observed in the Perseus and the Virgo clusters (Fabian et al. 2003; Forman et al. 
2005). Fabian et al. (2003) and Fabian et al. (2005) argued that the viscous dissipation of the 
sound waves is responsible for the heating of a cool core. They estimated the dissipation rate 
assuming that the waves are linear. However, when the amplitude of sound waves is large, 
the waves rapidly evolve into non-linear weak shocks (Stein & Schwartz 1972), and their 
dissipation can be faster than the viscous dissipation of linear waves 1 . Although Fabian 
et al. (2003) argued the presence of weak shocks, their evolution from sound waves was not 
considered. Ruszkowski, Briiggen, & Begelman (2004) performed numerical simulations of 
dissipation of sound waves created by AGN activities. Their results actually showed that the 
sound waves became weak shocks. However, their simulations were finished before radiative 
cooling became effective. Thus, the long-term balance between heating and cooling is still 
unknown. In this letter, we consider the evolution of sound waves to weak shocks, and 
analytically estimate the 'time-averaged' energy flux of the propagating waves as a function 
of distance, explicitly taking account of the dissipation at weak shock fronts and its global 
balance with radiative cooling. We assume the Hubble constant of H = 70 km s^ 1 Mpc -1 . 



1 Of course viscosity plays an important role for formation and dissipation of shocks. For simplicity 
however, we refer to the energy dissipation of linear sound waves through viscosity as 'viscous dissipation', 
and refer to that through formation of non-linear weak shocks as 'shock dissipation'. 
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2. Models 



We assume that sound waves are created by central AGN activities. The waves prop- 
agate in the ICM outwards. These waves, having a relatively large but finite amplitude, 
eventually form shocks to shape sawtooth waves (Landau & Lifshitz 1959; Mihalas & Miha- 
las 1984). If the velocity amplitude is larger than ~ 0.1 sound velocity (the Mach number 
is > 1.1), those waves steepen and become weak shocks after propagating less than a few 
wavelengths (e.g. Suzuki 2002). These shock waves directly heat the surrounding ICM by 
dissipating their wave energy. We adopt a heating model for the solar corona based on a 
weak shock theory (Suzuki 2002; Stein & Schwartz 1972, see also Fujita, Suzuki, & Wada 
2004). We assume that a cluster is spherically symmetric and steady. 

The equation of continuity is 



M = -Anr 2 pv , (1) 

where M is the mass accretion rate, r is the distance from the cluster center, p is the ICM 
density, and v is the ICM velocity. The equation of momentum conservation is 

dv = _GM(r) _ 1 dp 1 1 d 

dr~ r 2 pdr pc s {l + [(7 + l)/2]a w } r 2 dr 1 w> V> 

where G is the gravitational constant, M(r) is the mass within radius r, p is the ICM 
pressure, c s is the sound velocity, 7(= 5/3) is the adiabatic constant, and a w is the wave 
velocity amplitude normalized by the ambient sound velocity (a w = 5v w /c s ). The wave 
energy flux, F w , is given by 

F w = ^pcW w + • (3) 

The energy equation is 

d/1 j 7 k B T\ GM(r) 1 d 2 2 



p v T o v + — ~{ — + p° — 2^ + irn^ + F c)} + <^T) = , (4) 

dr \2 7 — 1 pm H J r l r l dr 

where ks is the Boltzmann constant, T is the ICM temperature, /i(= 0.6) is the mean 
molecular weight, inn is the hydrogen mass, n e is the electron number density, and A is the 
cooling function. The term V • F w indicates the heating by the dissipation of the waves. We 
adopt the classical form of the conductive flux for ionized gas, 

F c = -UoT 5 ^ (5) 
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with k = 5x 10~ 7 in cgs units. The factor f c is the ratio of actual thermal conductivity 
to the classical Spitzer conductivity. We adopt the cooling function for the metalicity of 
Z = 0.5 Z Q , which is the typical value in the central region of a cluster, 



n 2 e A(T) = [Ci(k B T) a + C 2 {k B Tf + C 3 ]mn. 



(6) 



where n,i is the ion number density, and the units for k B T are keV. The constants in equation 
(6) are a = -0.447, (3 = -0.232, d = 0.947, C 2 = -1.71, and C 3 = 0.922, and we can 
approximate riin e = 0.70(p/m H ) 2 . The units of A are 10~ 22 ergs cm 3 . This approximation 
reproduces the cooling function calculated by Sutherland & Dopita (1993) within ~ 10 % 
for 0.03 < T < 25 keV. The equation for the time-averaged amplitude of the shock waves is 
given by 



dr 



Oi w 

T 



1 dp 

p dr 



2(7 + IK 

C.,T 



1 dc s 
c, dr 
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where r is the period of waves, which we assume to be constant (Stein & Schwartz 1972; 
Suzuki 2002). The second term of the right side of equation (7) denotes dissipation at each 
shock front. We assume that the mass distribution of a cluster can be represented by the 
NFW profile (Navarro, Frenk, & White 1997), although recent studies suggest a little steeper 
profiles (e.g. Fukushige & Makino 1997): 



M(r) oc 



In 1 + 



r s (l + r/r s ) 



(8) 



where r s is the characteristic radius of the cluster. 



3. Results 

For parameters of our model cluster, we adopt the observational data of the Perseus 
cluster (Ettori, De Grandi, & Molendi 2002). We assume that r s = 280 kpc, M(r 1000 ) = 
3.39 x 10 14 M , and riooo = 826 kpc, where the mean density within ta is A times the 
critical density of the Universe. Waves are injected at the inner boundary r = r , which 
should be close to the size of bubbles observed at cluster centers. We assume that Ao = r , 
where Ao is the initial wavelength. If the waves are injected in a form of sound waves with 
amplitude 0.1 < a w < 1, waves travel about Ao before they become shock waves (Suzuki 
2002). Therefore, for r < r < r + A = 2A , we assume that V • F w = (eqs. [2] and [4]), 
and that the second term of the right-hand side of equation (7) is zero. The temperature, 
electron density, and wave amplitude at r = r are T n e0 , and a w o, respectively. Unless 
otherwise mentioned, the first two are fixed at T = 3 keV and n e0 = 0.08 cm -3 , respectively, 
based on the observational results of the Perseus cluster (Sanders et al. 2004). 
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In Figure 1, we show the results when r is fixed at 1 x 10 7 yr, and a w0 is changed. For 
the Perseus cluster, Fabian et al. (2003) estimated that a W Q ~ 0.5. The dissipation length 
is defined as l w = \F W /V • F w \. For these parameters, the initial wavelength is A = 9 kpc, 
which is roughly consistent with the Chandra observations (Fabian et al. 2003). The wave 
energy injection rate is given by ~ 47rrQ.F ?i ,(ro) ~ 10 45 erg s -1 , which is comparable to the 
jet power of the nucleus in the Perseus cluster (Fabian et al. 2002). Other parameters are 
M = 50 M Q yr" 1 , and f c = 0. In general, larger M reproduces observed temperature and 
density profiles better. However, large M is inconsistent with recent X-ray observations as 
was mentioned in § 1. For comparison, we show the results of a genuine cooling flow model 
(M = 5OOM yr _1 , a w0 = 0, and f c = 0) and the Chandra observations of the Perseus cluster 
(Sanders et al. 2004). Figures la and lb show that only a small region is heated. The jumps 
of T and n e at r = 2Ao = 18 kpc are produced by weak shock waves that start to dissipate 
there. The energy of the sound waves rapidly dissipates at the shocks, which is clearly 
illustrated in short dissipation lengths, l w ~ 2-15 kpc (Fig. lc). These dissipation lengths 
are smaller than those of viscous dissipation for linear waves, which can be represented by 
/„ = 420 Ag n .08 ^3~ 2 kpc, where the wavelength A = 9A 9 kpc, the density n = 0.08n 00 8cm~ 3 , 
and the temperature T = 3T 3 keV (Fabian et al. 2003; Landau & Lifshitz 1959). In Figure 1, 
the ICM density becomes large and the temperature becomes small at r > 2A so that the 
rapid shock dissipation is balanced with radiative cooling. Because of this, waves cannot 
reproduce the observed temperature and density profiles that gradually change on a scale 
of ~ 100 kpc. Note that the density peaks in Figures lb indicates that the solutions are 
unstable. In an actual cluster, this would lead to convection, and the ICM would be heated 
through the convection rather than the dissipation of sound waves alone. 

In Figure 2, we present the results when r = 2 x 10 7 yr. Compared with the case 
of r = 1 x 10 7 yr, the wave energy dissipates in outer regions. However, the dissipation 
lengths are still smaller than the cluster core size (~ 100 kpc). Note that larger r (or Ao) 
means formation of larger bubbles. As indicated by Churazov et al. (2000), it is unlikely that 
the size of the bubbles becomes much larger than 20 kpc; the bubbles start rising through 
buoyancy before they become larger. On the other hand, when r < 10 7 yr (A < 9 kpc), the 
waves heat only the ICM around the cluster center. The predicted temperature and density 
profiles are obviously inconsistent with the observations. 

The inclusion of thermal conduction changes the situation dramatically. Figure 3 shows 
the results when f c = 0.2. The value of f c is based on the study of Narayan & Medvedev 
(2001). The models including both wave heating and thermal conduction can well reproduce 
the observed temperature and density profiles. Figure 3c shows the contribution of the wave 
heating (—V • F w ) to compensating radiative cooling (rig A). Since —V • F w /n 2 e K > 1/2 
for r ~ 20-30 kpc, the wave heating is more effective than the thermal conduction in that 
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region. 

4. Discussion 

We showed that sound waves created by the central AGN alone cannot reproduce the 
observed temperature and density profiles of a cluster, because the dissipation length of the 
waves is much smaller than the size of a cluster core and the waves cannot heat the whole 
core. The same problem has been known for models of solar-corona heating by sound waves 
(e.g. Stein & Schwartz 1972). The problem of the short dissipation length should also be 
studied for other heating models including the bubble motion in clusters. On the other hand, 
we found that if we include thermal conduction from the hot outer layer of a cluster with 
the conductivity of 20% of the Spitzer value, the observed temperature and density profiles 
can be reproduced. The idea of the "double heating" (AGN plus thermal conduction) was 
proposed by Ruszkowski & Begelman (2002). 

However, the fine structures observed in cluster cores may show that the actual conduc- 
tivity is much smaller than that we assumed (e.g. Fujita et al. 2002; Mazzotta et al. 2003); 
the structures would soon be erased, if the conductivity is that large. If the conductivity is 
small, we need to consider other possibilities. While we considered successive minor AGN 
activities, some authors consider that rare major AGN activities should be responsible for 
heating of cool cores (Soker et al. 2001; Kaiser & Binney 2003). In this scenario, powerful 
bursts of the central AGN excite strong shocks and heat the surrounding gas in the inner 
region of a cluster on a timescale of > 10 9 yr. In fact, McNamara et al. (2005) found such a 
violent activity in a distant cluster. Moreover, in this scenario, heating and cooling are not 
necessarily balanced at a given time, although they must be balanced on a very long-term 
average. This is consistent with the fact that there is no correlation between the masses 
of black holes in the central AGNs and the X-ray luminosities of the central regions of the 
clusters (Fujita & Reiprich 2004). Alternative idea is that cluster mergers are responsible 
for heating of cool cores (Fujita, Matsumoto, & Wada 2004; Fujita et al. 2005). In this 
"tsunami" model, bulk gas motions excited by cluster mergers produce turbulence in and 
around a core, because the cooling heavy core cannot be moved by the bulk gas motions, 
and the resultant relative gas motion between the core and the surrounding gas induces hy- 
drodynamic instabilities. The core is heated by turbulent diffusion from the hot outer region 
of the cluster. Since the turbulence is produced and the heating is effective only when the 
core is cooling and dense, fine-tuning of balance between cooling and heating is alleviated 
for this model. 
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Fig. 1. — (a) Temperature, (b) density, and (c) dissipation length profiles for a w o = 0.3 
(dotted lines), 0.5 (thin solid lines), and 0.7 (dashed lines). Other parameters are r = 
1 x 10 7 yr, M = 50 M Q yr" 1 , and f c = 0. Filled circles are Chandra observations of the 
Perseus cluster (Sanders et al. 2004). The bold solid lines correspond to a genuine cooling 
flow model of M = 500 M Q yr^ 1 . 
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Fig. 2. — Same as Fig. 1 but for r = 1 x 10 7 yr (solid lines), and 2 x 10 7 yr (dashed lines). 
Other parameters are a W Q = 0.5, M = 50 M yr -1 , and f c = 0. 
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Fig. 3. — (a) Temperature, (b) density, and (c) dissipation strength profiles for a w o = 
and T = 3 keV (dotted lines), a w0 = 0.5 and T = 3 keV (solid lines), and a w0 = 0.5 and 
T = 2.5 keV (dashed lines). Other parameters are r = 1 x 10 7 yr, M = 50 M yr" 1 , and 
f c = 0.2. Filled circles are Chandra observations of the Perseus cluster (Sanders et al. 2004). 



